Decoherence of number states in phase-sensitive reservoirs 
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The non-unitary evolution of initial number states in general Gaussian environments is solved analytically. Decoherence in 
the channels is quantified by determining explicitly the purity of the state at any time. The influence of the squeezing of the bath 
on decoherence is discussed. The behavior of coherent superpositions of number states is addressed as well. 
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1 Introduction 

Recent developments in experimental cavity Quantum 
Electrodynamics and controlled atom-photon interactions 
allow for a direct investigation of deeply quantum me- 
chanical configurations of the light field. In particular, 
the deterministic production of low-order number states 
seem to be at hand, by means of micromaser techniques 
in high-Q cavities 1 1 1 and of strong coupling with trapped 
atoms |2|. The probabilistic generation of number states 
via conditional measurements and post-selection ||3l has 
been demonstrated as well, by exploiting parametric down 
conversion and low multiplication noise detectors |4|. A 
further possibility to generate number states with high fi- 
delities by atom-field interactions in high-Q cavities has 
been recently suggested l^. We also mention that, as 
for motional degrees of freedom, effective techniques to 
create number states have been developed and mastered 
16J. Moreover, the upcoming VLPC ('visible light photon 
counter') technology holds promising perspectives about 
the actual possibility of selectively detect low-order num- 
ber states |7 1. 

Besides being probes of fundamental quantum mechan- 
ical features, Fock states of the electromagnetic field are 
needed in several quantum information applications when- 
ever, for instance, reliable single-photon pulses are re- 
quired "91. These possibilities bring to the attention 
the problem of preserving the quantum mechanical prop- 
erties of number states, which are unavoidably corrupted 
by environmental decoherence. Indeed, their inherent non- 
classical nature makes such states especially fragile and dif- 
ficult to maintain. More specifically, the numerical analy- 
sis strongly suggests that the very possibility of generating 
pure number states is seriously hurdled by environmental 
decoherence, even in high-Q cavity settings IIIOI . In this 
paper we study the rate of decoherence of initial number 
states in general Gaussian noisy channels, i.e. in presence 



of dissipation in general Gaussian environments. The dy- 
namical behavior of the system is described by the quan- 
tum optical master equation, allowing for arbitrary phase- 
sensitive ('squeezed') baths. The dynamics will be solved 
in terms of the symmetrically ordered characteristic func- 
tion, while decoherence will be quantified by computing 
the purity fi — Tt during the evolution of the state. 

2 Solving the master equation 

Let us consider a denumerable, infinite-dimensional Hilbert 
space H, spanned by a Fock basis with n e N, of 

eigenstates of the hermitian operator h — a''a. The annihi- 
lation and creation operators a and at satisfy the canonical 
commutation relation [a, a^] — 1. We define the quadra- 
ture operators x = {a + a^)/^/2 andp = —i{a — a^)/V2 
describing, for instance, amplitude and phase quadratures 
of a single mode of the electromagnetic field, or position 
and momentum operators of a material harmonic oscillator 
Any quantum state of this system can be described ei- 
ther by its density matrix g or by its symmetrically ordered 
characteristic function x(a) 11 Tl . defined as 



X{a) = Tr{gDa) , 



(1) 



where Da = exp(aa 



t _ 



is the unitary displacement 



operator In the following we will make use of phase space 
variables x and p, defined by a = {x + ip) / \/2. Moreover, 
it is useful to define the covariance matrix cr, associated to 
a state g by 

<T,J^-{X,X,+X,X,)~{X,){X,), 

with xi = x,X2 = p and (O) = Tr {gO) for the operator 
O. The dynamics we will study can be modeled by the 
coupling with a continuum of oscillators, described by the 
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following interaction Hamiltonian 

( [W{w)a'^h{oj) + W{u))*ab'^] , (2) 



where stands for the annihilation operator of the bath 
mode labeled by the variable lo, whereas W{ijj) represents 
the coupling. The state of the bath has been assumed to 
be stationary. Under the Markovian approximation, such a 
coupling gives rise to a time evolution ruled by the follow- 
ing master equation (in the interaction picture) II12I 

Q = '1[n L[a'^]Q + {N + 1) L[a]Q - 

M* D[a]Q + M D[a)]Q\ , (3) 



where the dot stands for time-derivative, the Lindblad su- 
peroperators are defined as L[0]q = 20gO^ — O'^Og — 
QO^OwdD[0]Q EE 20gO-OOe-gOO,the coupling is 
7 — 2-kW'^{Q), while the coefficients N and M are defined 
in terms of the correlation functions {h^{Q)h{Lo)) — N5{lo) 
and (6(0)6(ijj)) — AIS{lu), where averages are computed 
over the state of the bath. The requirement of positiv- 
ity of the density matrix imposes the constraint |A/p < 
iV(iV + 1). At thermal equilibrium, i.e. for M = 0, N co- 
incides with the average number of thermal photons in the 
bath. If M Q then the bath is said to be 'squeezed', or 
phase-sensitive, entailing reduced fluctuations in one field 
quadrature. A squeezed reservoir may be modeled as the 
interaction with a bath of oscillators excited in squeezed 
thermal states |13|; several effective realization of such 
reservoirs have been proposed in recent years [14]. 

In general, the real parameter N and the complex pa- 
rameter M allow for the description of the most general 
single-mode Gaussian reservoir, fully characterized by its 
covariance matrix (Too, given by 



N+ RcM 
ImM 



ImM 

iV + Rc M 



(4) 



Such a Gaussian state constitutes the asymptotic state in the 
channel, irrespective of the initial condition. The param- 
eters 7, N and M completely characterize the Gaussian 
channel. A more suitable parametrization of the environ- 
mental state is provided by the following equations flSl 
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^J'c 



4IMJ 



cosh(2r) = ^n + 44jUM^' 



(5) 



(6) 



tan(2(p) = -tan(ArgMj) . (7) 

The quantities ii^o, r and Lp are, respectively, the purity, 
the squeezing parameter and the squeezing angle of the 



squeezed thermal state of the bath. This parametrization 
will prove useful in the following. 

Eq. (|3} is equivalent to the following diffusion equation 
for the characteristic function x in terms of the quadrature 
variables x and p 



X = 



X 

oo 1 __]X 



(8) 



It is easy to verify that, for any initial condition x{x, p), the 
following expression solves Eq. (|8} 

X = x{xe-i\pe-i')e-i^-P^--(iy'"^''"'^ . (9) 

The initial condition we will deal with is a number state 
|n)(n|, whose symmetric characteristic function Xn can be 
easily determined fl lj . and reads 



Xn{a) = {n\Da\n) = e L„{\a\'^) , 
where L„ is the Laguerre polynomial of order n 



Ln{x) 



m=0 



m 



(10) 



(11) 



Putting Eqs. (|9jl and ilO\ together and switching again to 
quadrature variables yields the solution Xn{t), accounting 
for the evolution of the initial number state | n) in the noisy 
channel 



Xn{t) 



with 



rrit) 



1 



-ft 



-It 



(12) 



(13) 



Notice that, clearly, a{t) is not the covariance matrix of 
the evolving state, because of the presence of the Laguerre 
polynomial in Eq. (I12> . 

3 Decoherence of number states 

Decoherence of the initial pure state in the channel will be 
quantified by following the evolution of the purity ji — 
Tr Q^. Such a quantity properly describes the degree of 
mixedness of a quantum state g. For continuous-variable 
(CV) systems it takes the value 1 for pure states (repre- 
sented by normalized projectors) and the value zero for 
maximally mixed states. The conjugate of /i is referred to 
as the 'linear entropy' Sl in information theory: Sl = 

The purity of a quantum state of a single-mode CV sys- 
tem is easily computed as an integral over the whole phase 
space L15J 



1 

2^ 



IxP da;dp 



(14) 
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The generalization of Eq. il4i to multi-mode systems is 
straightforward, and allows to track the dynamics in noisy 
channels of entangled two-mode Gaussian states 1 16| and 
of Schrodinger cat-like states 1 17|. Moreover, for CV sys- 
tems, it has been recently proved that knowledge of the 
global and marginal purities provides a strong and exper- 
imentally reliable characterization of the entanglement of 
arbitrary two-mode Gaussian mixed states |18|. In the 
present instance, Eqs. il2\ and (I14> allow to compute the 
purity fj,n{t) of an initially pure number state of order n 
evolving in the Gaussian channel. 

3.1 Thermal bath 

The instance of a reservoir at thermal equilibrium corre- 
sponds to the choice M = 0. In this case, both the envi- 
ronmental Gaussian state and the initial condition are rota- 
tionally invariant in phase space, see Eqs. and ilQ\ . so 
that polar coordinates constitutes a suitable and convenient 
choice. Employing the variable s = e~'''*|Q;p, one gets 

/•OO 

/,„(i)=e^*/ L2(s)e-(^'"(2^+i)-2^>ds. (15) 



Which can be analytically solved, using the relation 119\ 



[a - 2Y 



e--Li{x)dx='——^PAl 



a2 - 2a 



where P„ is the Legendre polynomial of order n 

1 d" 



Pn{x) 



2"n! dx" 



{x^ - IT , 



to find 



with ^ = e<\2N +1 



2N. 



(16) 



(17) 



Eq. ( II6I 1 provides the exact evolution of the purity of an ini- 
tial number state |n) in a thermal channel, fully determined 
by its mean photon number TV. Quite clearly, the purity in 
the channel is a decreasing function of N. It is a decreasing 
function of n as well: higher order number states are more 
fragile and decohere faster Moreover, number states with 
n > can show a local minimum of the purity and a partial 
revival up to the asymptotic purity — l/(2iV + 1). 

3.2 Squeezed bath 

We will now deal with the general instance M ^ 0. Due 
to the rotational symmetry of the characteristic functions 
of number states, we are free to choose any direction of 
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Figure 1: Evolution of purity for initial number states in Gaussian 
channels with /loo = 0.5 (purity revivals are evident in the plot). 
The solid line refers to state |l) in a non squeezed bath; the dotted 
line refers to state |l) in a bath with r = 1; the dashed line refers 
to state |2) in a non squeezed bath; the dot-dashed line refers to 
state |2) in a bath with r = 1. 



squeezing of the environmental state of the channel. In- 
deed, the rate of decoherence can only depend on the mod- 
ule I Af I of the parameter A/, which can therefore be chosen 
real and positive (corresponding to the choice (p = 0), with- 
out loss of generaUty. Making such a choice and exploiting 
Li9J 



2tt 



ef-^^d(^ = 2^/o(bl) 



(where Iq{x) = Jo{ix) = Y^kLo JWkTj^ ^he zero order 
modified Bessel function of the first kind), one eventually 
finds 



.7* 



-«^L„(s)/o(|M|(, 



,7* 



l)s) ds. 



(18) 

Such an integral cannot be further simplified, but can be nu- 
merically estimated to analyze the effect of the squeezing of 
the bath on the decoherence of number states. Eq. ( I18> ob- 
viously reduces to Eq. (I16> for M = 0. Besides, the asymp- 
totic expressions can be analytically integrated to yield 



lim finit) = Mo 



1 



v/(siV+ 1)2 -4|M|- 



which, according to Eq. (|5|i, is the asymptotic purity in the 
channel, fixed by the reservoir state, irrespective of the cho- 
sen initial condition. 

Eq. ( list shows that is an increasing function of \M\. 
However, the dependence on squeezing has to be analyzed 
by considering the parameters /ioo and r instead of N and 
AI, because they permit to study the effect of squeezing 
(quantified by r) at given asymptotic purity fioo- Such a de- 
pendence can be reconstructed inserting Eqs. (|5l|6l[n} into 
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Eq. ( list and turns out to be quite involved. Anyway, a nu- 
merical analysis has been carried out, and is summarized in 
Fig.^ where the evolution at short times (t < 7) is consid- 
ered. Note that this is the interesting time range, in which 
decoherence takes place before the system is driven to- 
wards the squeezed thermal state of the environment. Such 
an analysis definitely shows that squeezing the bath does 
not slow down the decoherence rate of number states in 
Gaussian channels. At a given asymptotic purity ^00 the 
highest purity is maintained for r — 0. 

4 Coherent superpositions 

We now consider the behavior of initial coherent superposi- 
tions of number states evolving in a general Gaussian noisy 
channel. To properly exempUfy the decoherence of such 
states, we focus on the simplest coherent normalized su- 
perposition iV'oi) = (|0) + e"' 1 1))/V2 (which constitutes 
a 'microscopic Schrodinger cat'). The characteristic func- 
tion xoi of this state is simply found lllTI 



2 (a e 



-iiS 



a e 



Xoi(a) 

(19) 

Inserting xoi as the initial condition in Eq. (|9} and perform- 
ing the integration of Eq. il4l yields, for the purity of the 
initial cat-like state evolving in the channel 

/ioi(i,r) - 4i.- e-27tZL(^^ + (eT*-l)(cosh(2r) 
+ cos(2i9 - 2(p)sinh(2r)) 

5 

+ e-^''*-^ (Voo + 8( e^* - l)/ioo cosh(2r) 
+ (e'^* - 1)2(3 cosh(4r) + 1)) (20) 



where 



-7* 



) +e 



-27* 



Mo 



■ cosh(2r) 



-1/2 



(21) 

is the purity of an initial vacuum in the channel 1151 . 
Eq. i2Q\ shows that, quite interestingly, the evolution of 
the coherent superposition is sensitive to the phase ip of 
the bath. It is straightforward to see that the optimal 
choice maximizing purity at any given time is provided by 
1} = (p + Tr/2. Fixing such a choice, we have numerically 
analyzed the dependence of /igi on r: for small squeezing 
parameters, the purity /ioi does increase with r. The opti- 
mal choice for r depends on time, for 'yt — 0.5 it turns out 
to be r ~ 0.28. The relative increase in purity is plotted 
in Fig. |2] as a function of time, for various choices of the 
squeezing parameter 7'. 




Figure 2; The relative increase in purity, defined by A/i/fi = 
{^ioi{t,r) - /ioi(i,0))//ioi(t, 0), as a function of time during the 
evolution of the superposition \tpQi) in Gaussian channels. The op- 
timal condition i9 = ip + n/2 is always assumed, while ^loo = 0.5. 
The solid line refers to a bath with r = 0.28, close to the optimal 
value; the dotted line refers to a bath with r = 0.4 and the dot- 
dashed line refers to a bath with r = 0.1. 



5 Comments and Conclusions 

We have analytically solved the non-unitary evolution of 
number states in Gaussian noisy channels, and quantita- 
tively estimated their decoherence. The dissipative model 
we have considered covers a variety of physical situations, 
from thermal dissipation of stationary modes in optical cav- 
ities to corruption of traveling waves in lossy fibers. In par- 
ticular, we have straightforwardly shown that Fock states 
of higher order decohere faster and that squeezed baths do 
not help to preserve the quantum coherence of initially pure 
number states. On the other hand, when considering coher- 
ent superpositions of number states, we have shown that 
squeezed reservoirs can help to slow down decoherence, 
provided that the phase of the bath is optimally locked to 
the coherent phase of the superposition, and that the inten- 
sity of the squeezing is properly chosen. This suggests, 
looking towards practical implementations, that feedback 
schemes (simulating squeezed reservoirs by quantum non 
demolition measurements) could indeed be helpful in pro- 
longing the lifetime of coherent superpositions of Fock 
states, i.e. number cat-like states. Moreover, it is worth 
noticing that the same effect on the preservation of the pu- 
rity of the superposition, warranted by squeezing the bath, 
can be obtained by an opposite squeezing of the initial state. 
In thermal baths, a superposition of squeezed number states 
proves to be more robust against decoherence than a super- 
position of non-squeezed number states. 
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